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Abstract 



It is widely believed that the pole mass of a quark is infrared-finite and 
gauge-independent to all orders in perturbation theory. This seems not to 
have been proved in the literature. A proof is provided here. 
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I. INTRODUCTION 



Long ago Tarrach showed that the pole mass of a quark is infrared-finite and gauge- 
independent through two-loop order in perturbation theory. It is widely believed that this 
result holds through any finite order in perturbation theory. There does not seem to be a 
reference in the literature containing a proof, however, and this paper aims to fill that gap. 
To be specific I consider QCD with one massive quark and Uf — l massless quarks, but extra 
massive quarks do not change the argument. 

There are a few reasons why one might suspect infrared divergences to arise in the 
perturbative-QCD series for the pole mass. QCD contains massless self-interacting gluons, 
so the infrared behavior is often worse than for QED. In QCD and QED, infrared divergences 
do arise in the two-loop self energy. Another worry is that the nonperturbative pole mass, 
if at all defined, is clearly very sensitive to the infrared. 

On the other hand, the kinematics of QCD, with the quark content under consideration, 
is like that of QED with a massive muon, some massless scalar bosons, and some massless 
electrons [0. Here (or with nf large enough to make QCD infrared-free) nonperturbative 
problems should not arise in the infrared. For the QED model one can even pick boundary 
conditions so that gauge-invariant muon states are in the spectrum.^ Also, the infrared 
divergence from the two-loop self energy is cancelled by a term from evaluating the one-loop 
self energy on the one-loop mass shell and iterating one could hope that this mechanism 
occurs at any order. Moreover, in asymptotically free QCD a remnant of the anticipated 
infrared sensitivity appears through renormalons 0,^; this implies that perturbation theory 
can "know about" the infrared behavior of nonperturbative QCD without having infrared 
divergences at fixed order. 

The above discussion is nothing but a duel of fears and hopes, and it should be replaced 
by definitive results. Below it is shown that the infrared divergences cancel in the pole mass, 
as at two loops, to any order. It is then a simpler matter to confirm that the pole mass is 
independent of the gauge chosen for the calculation. The arguments are straightforward and 
can be found in textbooks, although, to the best of my knowledge, the specific application 
is not. 

This paper is organized as follows: Some notation is in Sec. |I[ Section |T| proves that 
the pole mass is infrared-finite at every order in perturbative QCD. Section |V| proves that 
the pole mass does not depend on the gauge-fixing function through every finite order. Some 
concluding remarks are in Sec. 0. 

II. NOTATION 

This paper uses the metric g^u = diag(— 1, 1, 1, 1) in Minkowski space-time. In particular, 
the momentum p of a real particle with mass M satisfies = — M^. 

The pole mass is derived by identifying the pole in the massive quark's full propaga- 
tor S{p). One has 



The magnetic monopole considered by Ref. is a prototype of such a state. 
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(1) 



where mo is the bare mass and the self energy S is given by the sum of one-particle irreducible 
Feynman diagrams. One can write 

E(p, mo) = mo) + mo5(p^ mo), (2) 

exhibiting the parametric dependence of S on the bare mass mo- In considering the depen- 
dence on it is convenient to use a slight abuse of notation, S(p^,mo). 
The propagator S{p) has a pole at = — M^, where 

M = moZm (3) 

and 

Z„= Ihn '-^i";""; , (4) 

provided the limit is not infrared divergent. In perturbation theory one applies Eq. (^) by 
expanding the right-hand side through Lth order and setting to — iteratively. A gauge- 
invariant ultraviolet regulator, such as dimensional regularization or a lattice, is assumed 
but not made explicit. Because gauge theories are renormalizable, ultraviolet divergences of 
the coefficients are compensated by the bare mass and gauge coupling.^ 
Perturbative series are written, for example, as 

oo 

Zm = l + Y.9l'zt (5) 

1=1 

where is the bare gauge coupling. Below it is convenient to use a short-hand for 
the /th term in the perturbative series of expressions abbreviated here with •. 

The momentum p is reserved for the external momentum of the quark. Loop momenta 
are denoted generically by k. 



III. INFRARED FINITENESS 

To prove infrared-finiteness I follow the methods of Chapter 13 of Ref. [0. First, I 
recall how, in perturbation theory, one finds singularities in Green functions. This analysis 
establishes that the propagator and the self energy have a branch point at p^ = — M^. 
Since the pole mass requires the self-energy functions to be evaluated here, one must check 
whether they diverge at the branch point or not. It turns out that the on-shell self energy 
does suffer from infrared divergences, but I show that they drop out of the pole mass. 



^For example, in dimensional regularization one could introduce a renormalized mass m(fi) = 
moZjfi by minimal subtraction. Then one could focus on the infrared behavior of Zm/Z^. 
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A. Location of Singularities 



Consider an arbitrary Feynman diagram (of any Green function) , with quark propagators 
rationalized and all denominators combined with Feynman parameters a,. If the diagram 
has n lines the resulting denominator is 



D 



.1=1 



(6) 



where qi = qi{p, k) are the momenta of the internal lines. The Green function is an analytic 
function of Lorentz invariants of the external momenta, up to branch points. Branch points 
can arise only when D vanishes, but that is not enough. In addition, the contour of integra- 
tion (over Feynman parameters a and loop momenta k) must be pinched 0. This happens 
if and only if on each internal line 



2 I 2 



or aj = 



(7) 



(8) 



and, furthermore, following any closed path i in the diagram 

iee 

with the sign of q^ taken in the sense of the path i. Equations (|^ and (|) are the so-called 
Landau equations. 

Solutions of the Landau equations have a physical interpretation [Q. Up to an overall 
factor the Feynman parameter is the ratio of the time elapsed, from one end of the line 
to the other, to the energy propagating on the line. Thus, aiqi is the space-time separation 
between the two ends, and Eq. (P) says a loop in the diagram corresponds to a loop in space- 
time. Furthermore, Eq. (^ says that internal lines either are on shell [qf = —mj) or do not 
propagate (a^ = 0). For each diagram one obtains the reduced diagram by shrinking off-shell 
lines to a point. Then branch points arise if and only if the reduced diagram represents a 
genuine physical process of on-shell states. 

The physical picture given above is useful, because it is often easier to find solutions 
to Eqs. (|^) and (^ with physical reasoning instead of with algebra. For example, a two- 
point function has branch points only at normal thresholds, that is, when is just right to 
produce a collection of on-shell particles. For the massive quark propagator these branch 
points are at 



P 



'[{l + 2r)Mf 



(9) 



corresponding to creation of the massive quark plus r massive pairs. These branch points 
are accumulations of infinitely many solutions to the Landau equations, because once a 
solution is found, others are given by adding zero-momentum massless lines. Physically, this 
is because it costs nothing to create an extra soft gluon or extra soft pair.0 If the solutions 
accumulate too quickly, an infrared divergence will develop. 



^If one of the other quarks had a mass Mi, the branch points would be at 
2sM;]^, so infinitely many collapse to the same point as M/ 0. 



-[(l + 2r)M + 
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On the other hand, note that there are no colhnear divergences. As soon as the massive 
quark radiates non-zero momentum, it is off shell, and the (un) physical picture disallows a 
singularity. 



B. Infrared Divergences 

To examine the infrared properties, one performs a power-counting analysis. One scales 
some or all loop momenta by a factor A; if the Feynman integral scales as A'^ as A — 0, one 
says the degree of infrared divergence is /i. For example, in d dimensions the momentum- 
space volume element df'k has fi{d'^k) = d. If /i > 0, an integral is infrared convergent. 

The conclusions derived above for arbitrary diagrams apply equally well to the one- 
particle irreducible ones contributing to the self energy. It is convenient to route the external 
momentum p along the "main line," the massive quark line that runs all the way through a 
self-energy diagram. Off the main line the momenta are independent of p, and the degrees 
of infrared divergence are straightforward. Soft gluon and ghost propagators contribute 
/i(A) = —2, and soft massless quark propagators fJ,{S) = —1. Soft three-gluon and gluon- 
ghost vertices contribute //(Vs) = and other soft vertices fiiV) =0. In a closed loop the 
massive quark propagator So{k) = + mo) has degree of infrared divergence 0. 

The internal parts of the main line have propagator Sq{p + k). When k is soft 

' i{i)+^)+mo p^ + 2p-k + ml ^ ^ 

Off shell (away from the branch point) such lines have degree of infrared divergence 0. On 
shell, however ,0 

SoiP + k)^^^, (11) 

which gives degree —1. 

When all loop momenta are soft and p"^ = — M^, an arbitrary QCD (or QED!) self-energy 
diagram Gl with L loops, but no massive quark loops, has degree of infrared divergence 

fi{GL) = l + L{d-4) (12) 

in d dimensions. This holds at one loop. Higher-loop diagrams can be built by adding gluons, 
ghost loops, or (for now) massless quark loops. It is enough to insert the loops into gluon 
propagators; more gluons can be added later. Loop insertions give fi{GL+i) = fJ'iGi) — 4. 
New gluons give h{Gl+i) = fJ^iGi) + d — 2 + fii + /i2, where d comes from the new loop, —2 
from the gluon propagator, and the /ij are degrees associated with each end. The ends can be 
on a gluon or ghost line: /ij = 1 — 2 = — 1 for vertex and propagator; on a three-gluon vertex: 
/ij = — 1 from changing the three- to a four-gluon vertex; on the main line or a massless quark: 



^With = — one treats itiq — as higher order in Qq. 
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FIG. 1. Origin of the infrared divergence at two loops. 



Hi = —1 for the propagator.^ Thus, these ways all give fi{GL+i) = l^iGi) + d — A. Replacing 
an internal n- vertex polygon of a massless quark with a massive one increases ^{Gl) by n. 
Therefore, in four dimensions no infrared divergences can arise from the region with all loop 
momenta soft. 

Infrared divergences may come, however, from regions with some loop momenta soft and 
others not. According to the physical picture of the self energy, the non-soft lines can be 
shrunk to a point, augmenting the foregoing analysis with composite vertices. For n > 3, 
n-point vertices Vn contribute = to the total. Additional soft gluons attached to 

such vertices come with a propagator and a loop integration, adding d — 2io the total degree. 
Composite three-point vertices have the same infrared power counting as their fundamental 
counterparts. For example, gauge invariance guarantees the beneficial /^(Va) = +1 for the 
(composite) vertex of three soft gluons, in the same way it safeguards renormalizability. 
Thus, multi-point composite vertices do not lower the degree of infrared divergence in four 
dimensions. 

Internal, hard self-energy diagrams shrink to two-point vertices. In massless quark loops, 
the two-point vertex and extra propagator yield the harmless factor S(fc)/(i^) = A{k), as 
required by chiral symmetry. In massive quark loops the additional factor S(A;)/mo — > B{k) 
is also harmless. (These self energies are off shell and, therefore, well behaved.) Gauge 
symmetry provides two powers of soft momenta at two-point gluon and two-point ghost 
vertices, cancelling the extra propagator. Thus, these two-point vertices do not pose a 
problem. 

What remains are two-point vertices of the massive quark on the main line. The extra 
factor from inserting such a two-point vertex is S(p, mo)[mo — i'f\/2p ■ k, which lowers the 
degree of infrared divergences to /i < 0. For example, at two loops it is known |T|] that the 
self energy S[^l(— mQ,mo) is infrared divergent. The origin of the divergence, as the above 
analysis implies, is sketched in Fig. ^ The problem worsens at higher and higher orders, as 
more and more two-point vertices can arise on the main line. 



^When the massive quark is off sheU, the degree is increased by the number of main-hne propagator 
segments — even more convergent. 
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C. Infrared Cancellation 



For the pole mass at two loops the infrared divergence in {-ml, mo) is cancelled by 
the ©((^q) part of M^, mo). To examine this mechanism in general it is convenient to 

solve for the bare mass that implies a desired pole mass, namely 



mo 



MZ: 



with 



^-1 



1- A(-M2,M/Z„ 
l-5(-M2,M/Z„ 



(13) 



(14) 



showing the iterative nature of the solution. The Lth coefficient of the iterated 

expansion is infrared finite if and only if the coefficients Z]^, I < L, are infrared finite. 

The coefficients of the iterated expansion involve certain combinations of the coefficients 
of the self energy functions. Let 



s(p^Mz„^) = E^7o''sW(p^M) 



(15) 



define the coefficients ^^'^(p^, M) [and, by imphcation, A'^^^p'^, M) and B^^^p^, M)]. Iterative 
expansion yields 



l-l i 



1 



M) = M) + ^ ^ - [(Z-i - 1)^- 

i=i j=i J- 



dmi 



(16) 



mo=M 



Explicit calculation shows that T,^^^(p'^,M) = Yl^^^(p'^,M) is infrared finite and, when p'^ = 
— M^, gauge independent. At two loops S[^l(— M^, M) is also infrared finite and gauge 
independent, even though E[2](-M2,M) is not 
Since 



[z:, 



-11 [2] 



[1] 



(17) 
(18) 



with all self-energy functions evaluated at = — and mo = M, one has a basis for a 
proof by induction. Let us assume that the [Z-^Y^\ I < L, are infrared finite. To show that 
is also infrared finite it is enough to show that S[^l(— M^, M) is infrared finite. The 
induction hypothesis implies that S['l(— M^,M), / < L, is infrared finite (otherwise [Z^^]''' 
would not be). 

It would be a nightmare to identify all infrared divergences and verify cancellation on the 
right-hand side of Eq. (|TBp. Instead, it is more efficient to study S[^](p^,M) directly. The 
Dyson-Schwinger equation for the self-energy, depicted in Fig. ^ is a useful tool. Power- 
counting and the induction hypothesis together say that the only new infrared divergence 
at L loops can come from diagrams with the {L — l)th-order expansion of the quark propa- 
gator and an additional gluon: 
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FIG. 2. Dyson-Schwinger equation for the self energy. Grey blobs denote full propagators and 
white blobs denote one-particle irreducible functions. 



^IR 



(27r)4 P 



7;.[^(p + A;)][^-%. 



(19) 



Diagrams with a higher-order IPI vertex function or gluon propagator would have a quark 
propagator with L — 2 (or fewer) loops; they can be infrared divergent only if S^, for some 
/ < L, were too — contrary to the induction hypothesis. 

The key to obtaining from the right-hand side of Eq. ( [T9| ) is to write {q = p + k) 



Siq) 



li + M- [S(g, MZ„i) - M(Z^i - 1)] 



treating + M as 0((y'o) and expanding the bracket (iteratively) in g^. This expansion, for 



-M and k soft, produces a sum of chains 



[S{p+k)] 



X-l] 



M - 
2p-k 



Y^l[---\i:ip + k,M)-M{z^'-i] 



[ij] M -i^ 
2p-k 



(21) 



In any term of the sum the factors' superscripts Ij add up to L — 1, and the sum is over all 
such partitions of L — 1. Since k is soft the quantity in brackets reduces to 



E(p + k,M)-MiZ-'-l) 



V] 



[I] 



t^A{p\ M/Z^) + MZ~^B{p\ M/Z.^) - M{Z;^^ - 1 
ii) + M)A^^\-M'^,M) + 0{k). (22) 



The second step follows from Eq. ( |l4|) and setting — M^. In Eq. ( |2ll) the right-hand side 
of Eq. (p2D multiplies M — z^; the 0(1) part of the product vanishes for p^ = — M^, leaving 
a remainder of O(fc).0 Consequently, [S{p + k)]^^^^^ has degree of infrared divergence — 1, 
just like So{p + k). Thus, S[^](— M^, M) is infrared-finite, as was to be proved. 



^One might worry whether the remainder is an infrared-divergent derivative of the self energy. 
But because p -|- A; is off shell by /c, the remainder is proportional to A; In A;. 
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Although the above formulae are a bit clumsy, the mechanism that cancels the infrared 
divergences is simple. Equation (^OD says to split the bare mass into the pole mass plus 
a counterterm. The counterterm, like the shrunken self-energy subdiagram, produces a 
two-point vertex. Infrared divergences cancel in the (next order's) pole mass, because the 
combination of the two does not degrade the infrared power counting. In QED, this mech- 
anism was identified in a footnote to Ref. Il9|. 



IV. GAUGE INDEPENDENCE 

The gauge invariance of the mass renormalization factor is "nearly obvious." If the 
ultraviolet regulator respects gauge symmetry, the bare mass has a gauge-invariant meaning. 
From a physical point of view it would be unsettling if the pole mass were to depend on 
the gauge. Thus, the ratio = M/rriQ ought to be gauge invariant too. Without the 
infrared-finiteness established above, however, a proof would require painstaking separation 
of infrared and ultraviolet regulators. 

With infrared-finiteness one can study the gauge dependence by treating the massive 
quark like a normal particle, as long as one discusses perturbation theory only. For 
near — the propagator takes the form 

where Z2 is the field renormalization factor.[] Because the pole mass is only defined pertur- 
batively, M denotes here the pole mass through some finite order in perturbation theory. 

The previous section assumed a gauge fixing term (A/2) /(9 ■ A^d'^x. [In Eq. (p!9[), 
^ = 1 — A^^.] Suppose M depends on A. A shift A A induces a first-order change 



A 0/ N 

AS{p) 



dZo DM Zo 



dX dX li) + M 



(24) 



Note that dM/dX multiplies a double pole. 

On the other hand, the propagator is given by 

S{p) = FT{0\T^{x)ij{y)\0). (25) 

From Eq. (^), the change is 

AA r 

ASip) = —FT J d''z{0\Ti;ix)^P{y)id ■ Aiz) f\0) (26) 

A double pole would develop if the quark could scatter off (two) gluons with scalar polar- 
ization, which cannot happen because the scalar polarization decouples from the physical 



'^Z2 is not infrared divergent when is close to, but not equal to, — 
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state space. The mass M cannot, therefore, depend on the gauge parameter, although the 
residue Z2 certainly can. 

(Conversely, one can see immediately that the insertions generated by a shift in the bare 
mass or gauge coupling would develop a double pole and, thus, a shift in the pole mass.) 

For a more general gauge- fixing function f'^{A{z)) (and change A/°) the argument is 
similar. Let s denote the BRS operator and r] {ff) the (anti)-ghost field. The change in the 
two-point function involves [|T^] 

{s[Uz)^r{zmx),l,{y)) ~ {Uz)M\z)s[^{x)^{y)]) (27) 

~ {n,{z)Anz)t\r^,{x) - r^,{y)]^{x)i,{y)) (28) 

As before, but now because the ghosts decouple, this expression cannot develop a double 
pole. 

V. CONCLUSIONS 

The pole mass is widely used in the phenomenology of QCD and, when quark momenta 
are small compared to the mass, in nonrelativistic QCD and heavy-quark effective theory. 
In many of these contexts it is natural: it has considerable intuitive appeal, and it can be 
calculated with any ultraviolet regulator and in any effective theory. In some circles the 
pole mass — as an experimental quantity — has rightly fallen into disfavor, because infrared 
renormalons obstruct an unambiguous determination [^0. Had the pole mass turned out 
to be either infrared divergent or gauge dependent, one ought to have abandoned the pole 
mass for more basic reasons. 



Fortunately, Sec. |T| shows that the perturbative pole mass is infrared finite, even though 
the on-shell self energy is not. At higher orders the pole mass requires an iterative expansion 
and, thus, (infrared divergent) derivatives of the self energy. In the pole mass the total 
infrared divergence vanishes. The cancellation mechanism and, indeed, the power counting 
in QCD are the same as in standard QED. Attaching a virtual photon to an on-shell massive 
electron line has the same effect as attaching a gluon anywhere in a massive-quark self- 
energy diagram. I cannot imagine that infrared-finiteness of the electron mass has never 
been proved, but, except for a footnote §], I have not found a reference with a proof. 

Because of its physical appeal, the pole mass remains valuable theoretically. In addition 
to matching to effective theories, mentioned above, it is similarly useful in relating lattice 
QCD to continuum renormalization schemes [|1T|. A example of considerable phenomeno- 
logical interest is the application of (NR)QCD to threshold production of heavy quarks. 
There the pole mass is nearly irresistible, but it has been pointed out recently that infrared 



sensitivity in the mass's definition is conferred on the QCD potential as well |[T^-[14|. In- 



deed, Ref. [0 notes that a formula equivalent to Eq. (|16D would be helpful in showing that 



infrared divergences in the pole mass cancel at every order. 
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